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We formulate, with full generality, the asymptotic estimation theory for Gaussian states in terms of 
their first and second moments. By expressing the quantum Fisher information (QFI) and the elusive 
symmetric logarithmic derivative (SLD) in terms of the state's moments (and their derivatives) 
we are able to obtain the noncommutative extension of the well known expression for the Fisher 
information of a Gaussian probability distribution. Focusing on models with fixed first moments 
and identical Williamson 'diagonal' states -which include pure state models-, we obtain their SLD 
and QFI, and elucidate what features of the Wigner function are fundamentally accessible, and at 
what rates. In addition, we find the optimal homodyne detection scheme for all such models, and 
show that for pure state models they attain the fundamental limit. 



PACS numbers: 03.65.Ta, 03.67.-a, 06.20.Dk, 42.50.St 

Estimation theory plays a central role in modern devel- 
opments of quantum enhanced metrology. From a practi- 
cal point of view, it allows to assess the ultimate precision 
limits of given metrological schemes. From a fundamen- 
tal perspective, it provides a gold standard upon which 
to asses distinguishability of quantum states. Quantum 
estimation theory is an old subject (T], Q and it has seen 
huge developments over the last 20 years It has 

played a major role in understanding the fundamental 
powers and limitations of quantum measurement. The 
groundbreaking advances in atomic clocks and precision 
metrology 0| suggest that quantum estimation theory 
will only become more relevant as measurement precision 
reaches its fundamental limits. 

Common to almost all the disciplines of physics where 
precision metrology can provide significant results, is the 
fact that they benefit from the simplicity and power of 
the Gaussian state formalism Q. The latter has already 
proven its success and serves as an invaluable tool in 
describing quantum states of light and atomic ensembles, 
as well as providing useful insight and intuition. 

Despite the great success of both the Gaussian state 
formalism, and quantum estimation theory, these two 
have never been successfully merged. Indeed, the non- 
trivial equations defining central objects in quantum es- 
timation theory often forces to numerical methods [9l-[ll| 
for computing precision bounds and determining opti- 
mal measurements, and these difficulties are only aggra- 
vated by the infinite-dimensional nature of bosonic sys- 
tems. Most remarkably, these difficulties are not allevi- 
ated by the Gaussian state formalism, but for very par- 
ticular cases and without explicit harvest of the phase- 
space structure. Ironically, it is in the seminal book by 
Holevo where one encounters one of the first com- 
prehensive accounts of the phase-space formalism and 
Gaussian states on the one hand, and the foundations 
of quantum estimation theory on the other. There, cer- 
tain essential features of Gaussian states regarding opti- 
mal detection were established, but the analysis focused 
on the Gaussian shift model. This model, despite being 
extremely relevant, and certainly the first candidate to 
be studied, is not well suited for modern entanglement- 



enhanced metrology, where the signal is encoded in the 
state's correlations rather than the amplitudes. It is the 
purpose of the present work to provide a fully general 
phase-space formulation of the central quantities in quan- 
tum estimation theory, namely, the symmetric logarith- 
mic derivative (SLD), and the SLD quantum Fisher in- 
formation, with focus on general Gaussian states. After 
presenting this new formulation, we use it to study a wide 
class of models which include all pure state models. We 
also address the optimality of Gaussian measurements 
and derive sufficient conditions under which they are op- 
timal. 

Given a quantum model, i.e. a parametrized set of 
quantum states S = {pe}i the quantum Cramer- Rao 
bound (QCRB) establishes a lower bound to the vari- 
ance of any unbiased estimator of parameter 6, 



{Mf>iQ{eY 



(1) 



where Iq{0) is the SLD quantum Fisher information 
(QFI). The QFI is defined in terms of the symmetric 
logarithmic derivative (SLD), which is the Hermitian op- 
erator Cr that satisfies 



dope = + Co 



and the QFI reads 



(2) 



(3) 



Although the QCRB only establishes a bound, for uni- 
parametric models it is asymptotically attainable (l2| - 
[l4| . Therefore, it establishes the best asymptotic rate at 
which statistical fluctuations can decrease when measur- 
ing parameter 9 [J-Q. Thus, it is a quantitative measure 
of distinguishability of a state pg from its neighbors pe±s9 
and as such, is intimately related to the quantum Fidelity 
and the Bures distance [la]- On the other hand, the SLD 
not only has a geometric meaning; it also represents, by 
construction, an optimal observable [H, [l^l , in the sense 
that -to leading order- {Cg)g+se is proportional to the 
deviation from a reference state pe, 56, and it minimizes 
the statistical fluctuations. 
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We focus our analysis on systems of n bosonic 
modes, described by Hilbert space H — L'^{W'). 
These systems are characterized by 2n canonical vari- 
ables i?' = (Qi , . . . , Qn, Pi , . . . , Pn) with canonical com- 
mutation relations [K^jW] = iu^^ , where w is the 
2?! X 2n symplectic matrix, and il is the degenerate 
symplectic inner product, il{a, b) = '^ij o.^tf'^''^ = 
— aj3- The Weyl (displacement) operators are de- 
fined as WiO =expil7(^,i?), so that W''{^)K'W{C) = 
R'+C, and W{^)W{t]) = W{C + Tj)ei^^^''^\ In ad- 
dition, we introduce the symmetric product A o B — 
{AB + BA)/2. The symmetric product is not associative. 
We define it to have precedence over ordinary product, 
AoBC =iAo B)C ^Ao (BC). 

Gaussian states are defined as those states p e 
B^y.) having Gaussian characteristic function XpiO ~ 
tT[pW{^)]. Let A42n be the space of 2n x 2n matrices 
over M. We define the first and second moments as 

d' = tr[i?V] (4a) 
r*^' = 2tr[(P* - d') o {R^ - d^) p] , (4b) 

where d & B?'' and lo ,T e X2n(M). With these defini- 
tions, a Gaussian characteristic function reads 

xp{o^c^v{^ed-\en) (5) 

where d = —Lud and f — loT tu^ . 

The main goal of this work is to provide a formulation 
of Eqs. ([2]) and ([3]) in terms of dg, Fg, dgde and dgVo, and 
to explore the benefits of such formulation. Let S = {pg} 
be a Gaussian model with parameter 9. Any such model 
is fully described by the first and second moments dg 
and Tg of pg. As can be expected from the structure 
of Gaussian states, -and we show in Appendix [Xl - the 
SLD is quadratic in the canonical operators. It has also 
zero expectation as follows from its definition Eq. ([2]) and 
(?tr[pe] = 0. Taking as ansatz for Cg the expression 

ij 

+ Y.h{R -d^)-\r[LT] (6) 

i 

in Eq. © and taking the characteristic function of 
both sides one can relate the derivative of Xp(C) w.r.t 
to 9, dgXpiO (Ihs) to expectations of i?' o W{£,), and 
(i?* o R^) o W{£,) (rhs). This relation in turn, implies a 
relation between d, F, dd and dT (Ihs) and suitable lin- 
ear combinations of dijW{^) (rhs), where di = 9/i9^'. 
Equating terms with equal powers of ^ yields solutions 
to L and b in Eq. We report the solutions here and 



^ Summation of latin indices i and j runs from 1 to 2n, except for 
symbol k which labels the n modes and thus runs from 1 to n 



the technical details in Appendix [X] Define the linear 
map Vx ■■ M2n M2n as 

VxiY) ^ XYX'^ ~ ujYuj'^ . (7) 
With this, L G M2n and b G M^" are given by 

b = 2r-^dd, (8a) 
L = V^\dr), (8b) 

where the a pseudoinverse (Moore-Penrose inverse) is un- 
derstood whenever Vr is singular. The map 2?p^ in 
Eq. (j8bp will play a central role in our discussion. No- 
tice that VriX^) ~ VpiX)^ is symmetry preserving. 
Let Y = 2?p^(9F). Then Y satisfies the Stein equation, 
Y - FFF^ = -5r-i, where F = (Tuj)-^ and its spec- 
trum is contained in the unit disk (|_F| < 1). Assuming 
p is nonsingular, \F\ < 1 and the unique solution to the 
Stein equation is |16| 

oo 

Y = - ^F'^dT-^F^''. (9) 

A more detailed analysis of the structure of Vr and its 
pseudoinverse is given in Appendix |Al The following re- 
lation will be useful, 

dT = -Vr{dT-^) - ljOT-^oj'^ , (10) 

as follows from dT = -FaF-^F. 

As a first observation regarding the structure of the 
SLD obtained in Eq. ©, setting i9F = in Eq. ([Sb]) 
recovers the Gaussian shift model, already studied ex- 
tensively in the literature 0, [13] ■ In this case we obtain 
Cg = dd-T~^ ■ R, linear in the canonical operators, hence 
recovering the optimality of homodyne detection. 

Beyond the Gaussian shift model, for generic Gaussian 
models the SLD is at most quadratic in the canonical 
operators. Defining operators i?* = i?* — + i[L^^6]* 
one can write 

ij 

where C is a scalar. In addition, as follows from Eqs. (|8b|) 
and ([H]) the matrix L is the image of —dT~^ under the 
action of a completely-positive map. Therefore, when- 
ever i9F^^ or —dT~^ is positive semidefinite, L is so 
too, and there is a symplectic transformation Tg such 
that L = rg'^diag(Q!i, ai, ... , a„)Te. Thus, defin- 

ing i?* = J2j R^ and the number operators, Nk = 
\[{R^Y + (i?"+'^)2 - 1] = \[{Q''f + {P^f - 1] one has 

n 

Cg = Y,ak{Nk^ {Nk)g). (12) 
fc=l 

Therefore, in the case where dT~^ has definite signature, 
the SLD reduces to photon counting in suitably defined 
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modes. This result extends and generahzes the findings of 
[3 [3 J where the SLD was shown to have this structure 
for certain classes of channel estimation problems. Still, 
despite having a physical interpretation of the operator 
Cg (Gaussian transformations and photon counting) , im- 
plementing a measurement of it may still be prohibitive, 
especially if strong squeezing of the signal is required. 

Before considering more practical measurements, let 
us obtain the fundamental limit to their performance. 
From Eq. (O one can readily obtain the expression for the 
quantum Fisher information for general Gaussian mod- 
els. It is convenient at this point to endow M.2n with 
an inner product structure, (-^^1^) = tr[X^F], thus re- 
garding matrices as vectors (kets) \X), and linear maps 
thereof as B\X) = \AXB^), so that we can write 
2?r = F ® F - w ® w and Eq. (HbJ reads \L) = V^^\dT). 
With this inner product, I?r and its (pseudo)inverse are 
self-adjoint. The QFI reads [see Appendix IB] for details] 

Iq = i(o>F|(F®F - ijj®oj)-^\dV) + 2dd^T-^dd. (13) 

This expression is, for the first time, the most general 
form of the QFI in the Gaussian state formalism, and 
together with Eqs. ^ and ^ constitutes our main re- 
sult. It allows to compute precision bounds for a number 
of situations, and expresses the Fisher information in a 
form amenable for numerical computation, overcoming 
the difficulties posed by the infinite-dimensional char- 
acter of Eq. Notice that putting the dimensions 
back in and taking the classical limit {h 0) yields 
(T®T -h^ui®uj)-^ F"i(8)F~\ thus recovering 



2dd^T-^dd, (14) 



the Fisher information of a Gaussian probability distribu- 
tion centered at d with covariance Fg. Thus, Eq. (fT3|) is 
the noncommutative generalization of Eq. p4p . Indeed, 
^ plays an essential role in capturing the geometry and 
distinguishability properties of Gaussian states, and the 
term proportional to }t?uj®w accounts for the uncertainty 
due to noneommutativity of the canonical observables. 

Eq. ([9]) gives the unique solution for models with non- 
singular p, namely, no vacuum modes in the Williamson 
decomposition. This rules out the important case of pure 
state models. However, a generic solution to the Stein 
equation can be given for the class of models with co- 
variance matrix satisfying the relation 



(15) 



where i/ oc / is constant and we will treat it as a scalar. 
We call these models isothermal due to the constant 
temperature of their Williamson decomposition. These 



models include all pure state Gaussian models, as well 
as some mixed models that have recently attracted at- 
tention [i^l- It is easy to check that Eq. (fT5)) implies 
(i/2 - l)r = j/2l)p(F^i), which combined with Eq. pU)) 
leads to (\+v'^)V^^ [dV) = -v'^dV-^ . Despite some tech- 
nicalities in simplifying 2?p ^ oPr when Pr is singular, the 
resulting expression is also valid for pure state models, as 
can be shown by a detailed analysis of the singular case 
[v ^ 1). The QFI for such models is readily obtained 
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2 1 + 1/2 



tr[9FF-iaFF-i] + 2dd^T-^dd. (16) 



First, notice that the contribution due to first moments 
is equal to that of Eq. (|14p . This is due to the fact that 
there is always a reference frame for which the derivative 
dd is along a set of mutually compatible variables. In ad- 
dition, setting [v = 1) recovers the known expression for 
the QFI for pure state models [2l| . Interestingly, the cor- 
rection factor for thermal models v"^ + approaches 
1 in the large temperature limit (y ^ 1), recovering the 
Fisher information contained in the Wigner distribution. 
As noneommutativity of the canonical variables dictates, 
a faithful sampling of the Wigner distribution is beyond 
reach except in the high temperature regime, when ther- 
mal fluctuations render quantum fluctuations irrelevant. 
This explains why, in the high temperature limit, Iq ap- 
proaches the Fisher information of the Wigner distribu- 
tion. In addition, notice that for pure state models with 
fixed first moments {dd = 0), Iq is exactly 1/2 of the 
Fisher information contained in the state's Wigner func- 
tion. This fact deserves further attention. Consider the 
the model (F, dT) in the Williamson form 



Ffl — SrvSo 



(17) 



The essential quantity in Eq. (HH) is W = Sg^dYgSg'^^ , 
which is nothing but dVg expressed in the coordinates for 
which Tq is diagonal. However, the canonical transfor- 
mation Se is only determined up to an orthogonal sym- 
plectic transformation. Conveniently, W is symmetric 
and Hamiltonian [W o u = 0), thus there is always a 
symplectic orthogonal transformation O that diagonal- 
izes it. Therefore, defining the symplectic transforma- 
tion T — OS~^ and corresponding canonical coordinates 
iQi,...,Qn,Pi,...,Pn) = R = TR, one has 



TFT' = 



f 
ly 



TdTT' = 



where A > and we have made the v dependency ex- 
plicit in TdTT^ . This illustrates a characteristic trait of 
all models of the kind p5|) . i.e., that there exist canoni- 
cal coordinates W for which variations of the parameter 
correspond to a collection of single-mode squeezing oper- 
ations. In addition, it is clear that 



^ The factor of 2 in the second summand is due to the factor 2 in 
the definition of the covariance matrix, Eq. l l4bll 



(19) 
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Now consider the class of homodyne measurements ob- 
tained by measuring quadratures {Qk} (or {Pk})- The 
outcomes of such measurements are Gaussian distributed 
with covariance matrix = vl and d^g = ±v\, hence 
yielding /d given by 



(20) 



One immediately sees that for pure states [v = 1) this 
is optimal, as follows from Eq. (|19|1 . Iq = Id- Hence, 
no other measurement can perform better. Let us pause 
for a moment to discuss what we mean by homodyne de- 
tection in this general multimode setting. In principle, 
we regard homodyne detection as the measurement of 
any set of compatible canonical variables, labelled {Qk} 
for a suitably chosen canonical coordinates, R = TR. 
However, the set of variables {Qk} is only a convenient 
way to account for the information that the measure- 
ment outcomes provide. The real estimator O is some 
linear combination of the outcomes, and thus can be de- 
scribed as a linear combination of the canonical observ- 
ables = R, where a = (ai, . . . , a„, 0, . . . , 0) £ R^". 
Expressing O in terms of the original coordinates we 
get = R, with a ~ T^a. Hence, it remains to 
be seen that any such linear combination R can be 
implemented by simple linear combinations of the origi- 
nal Qk quadratures or passive Gaussian transformations 
thereof, and does not depend critically on the simulta- 
neous measurement of incompatible variables or on some 
impractical active transformations. Parametrizing pas- 
sive transformations as orthogonal symplectic matrices 
V e M2n; we seek V such that = (Va)'^ i^R), where 
a = Va ~ (.9,0). Writing V with blocks c,s £ Ain we 
get 



Va 



c s 

— s c 



car- 



(21) 



is al- 



so the condition that [X^a]p ~ — sc 
ways achievable by e.g., c = diag{aq-^k/\c(q:k + icep-k\) and 
s = dia.g{ap-k / \(Xq;k + 'ictp:k\)- Let g = caq + sap, then 
= a^ R = J2k9kQk [R = VR] shows that any linear 
combination of the canonical operators is implementable 
by passive transformations and homodyne detection on 
the original n modes. 

Going back to Eqs. pQ]) and ([20|l . it is clear that for 
mixed states (v > 1) there is, potentially, room for im- 
provement, as Iq > Id. More general Gaussian mea- 
surements can be implemented by attaching a Gaus- 
sian ancilla and performing homodyne detection [i^, [1^ . 



Adding an ancilla corresponds to replacing F — )• F © 7 
and dT 9F © 0, thus A A © 0. Therefore, with- 
out loss of generality we can focus on homodyne detec- 
tion. Consider the Fisher information Id of observables 
Q = {Qi, ■ ■ ■ ,Qn} out of a generic set of canonical co- 
ordinates i?* = J2j U"^^ RP , parametrized by a symplectic 
matrix U 



U = 



(22) 



Outcomes q = {gi, . . . ,q„} arc distributed according to 
the covariance matrix 7 = i'[UU^]qq ~ ^{aa^ + bb^). 



where the 



J 99 



subindex indicates that the block corre- 



sponding to Q quadratures has to be taken. Likewise, 
the derivative ^7 is of the form d"/ = v{a\aJ — bXb^). 
One can check that this model yields Fisher information 



Id{a,b) = itr 



(0,(A)-0b(A))' 



(23) 



where 4>zix) = 



-1/2 



T ^ — 1 /9 
ZXZ^ ' , 



with 



G {a,b}. The 

maps 4>a,b are completely positive and = 0a + is uni- 
tal. Using the property ab^ = ba^ one can show that </> is 
also trace-preserving. Finally, adding 2tr[0a(A)0b(A)] > 
to Eq. ^ yields Id{a,b) < itr[(/.(A)2], and us- 
ing Kadison-Schwartz inequality one gets Id{a,b) < 
^tr[0(A^)] — /*[ because (j) is trace-preserving. 

This shows that, for models of the form of Eq. (|15|1 . 
with fixed first moments, homodyne detection of a suit- 
able quadrature is always optimal among Gaussian mea- 
surements, with Fisher information upper bounded by 
/*[, and heterodyne detection cannot improve its perfor- 
mance. Moreover, for pure models it is optimal in a wider 
sense -among all quantum measurements-. The optimal- 
ity of homodyne detecti on g eneralizes and puts in context 
some earlier results [21I [2J] . 

Further work can be envisaged in various directions. 
On the more practical side, identifying the most general 
class of models for which homodyne detection is optimal. 
On the theoretical side, application of these results 
to the study of Gaussian channels is a natural way of 
proceeding. In addition, one expects that some extension 
of the methods used here may be able to provide insight 
to paradigmatic non-Gaussian models such as phase 
diffusion and degaussified states. 

We thank H. Cable, M. Hayashi and B.-G. Englert for 
useful discussions. The Centre for Quantum Technologies 
is funded by the Singapore Ministry of Education and the 
National Research Foundation as part of the Research 
Centres of Excellence programme. 
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Appendix A: The symmetric logarithmic derivative 

We will make use of Einstein's summation convention. 
In addition, we will make a distinction between covari- 
ant and contravariant indices to bookkeep the transfor- 
mation rules to which they comply to. In this spirit, the 

so 



inverse of the symplectic metric a; is flij = [ 
that QijUj^'' = 6f. We give a different symbol to avoid 
confusion, because componentwise i^ij ~ — w*^ . Also, we 
define di = and d with no index refers to d/dO. 

Taking the trace of Eq. ^ with the Weyl operators we 
get 

deXp{0=^^[i^BopW{0]. (Al) 
As an ansatz, suppose Co is at most quadratic in the 



canonical operators 

so the RHS of Eq. (|A1[) can be written as 

tT[jCg o p W^O] = L^°^ + ir^tr[p R' o WiO] 

+ L[ftr[p{R'oR^)oWiO] 



(A2) 



(A3) 



On the other hand, the Gaussian characteristic function 
allows to write 

deXpiO = {^idgd- \eden)xp{£.)- (A4) 
where d and T in proper covariant-contravariant notation 



read di = Qijd^ and Tij 



We begin by developping the RHS. The Weyl operators 
are defined as 



W{0=exp{in{C,R)), 

and we can show 
or equivalently 



(A5) 
(A6) 
(A7) 



It is tempting to define lower-indexed operartors Ri — 
riijR^ so that one can write diW{^) = iRi o W{£,). How- 
ever we will avoid this notation since the metric fl is an- 
tisymmetric and it would be essential to be consistent on 
which index R^ is contracted with. The second derivative 
reads 



^^^JW{0 = -^^knJl R'' O (i?' O WiO)- 



(A8) 



The symmetric product is not associative (A o B) o C ^ 
A o (B o C). To put Eq. (|A8P in a manifestly symmetric 
form we use 

R' o {R' o w{0) = (i?' o R') o w{o - \e^'w{o 

(A9) 

so that Eq. (jA8j) reads 

^^,w{0 = -n^k^ji (^R'' o R' - i^'^e') o W{0 (AlO) 

which is manifestly symmetric under From this we 
can derive the relation 

{R'' o i?') o W{0 = Qe''^' - u;'''cj'^9y^ iy(0- (All) 
Combining Eqs. (|T7l) and ([XTT]) with Eq. we get 

tT[po£.gW{0] (A12) 

= [L^'^-^L['^u;^'^^, + L^^\\ee-u^V^u))xp{^)■ 
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Using Eqs. (M^i and (|A12p . Eqs. ((XT|) reads 
(iC5(i-i^^af^)Xp(e) (A13) 

It is straightforward to evaluate the derivatives of Xp- 
These are 

(A14) 

and 

Lifilee - c.''=c.^'afeOxp(0 = (A15) 



Rewriting Eq. (jA13[) . and recalhng that Xp is nowhere 
zero, we get 



(A16) 



Equahng the different orders in ^ we get 

+ L(i)d + dTi(2)rf + itr[i(2)r] = (Al7a) 

L^^') + {L^^) + L^^)^)d^2r-^dd (AI7b) 

TL^^^r - ujL^^^uj'^ = 9r. (Al7c) 

By defining the hnear map Vx as 

VxiY) = XrX^ - wFw^, (AI8) 

Eq. (|A17cP can be written as 

2?r(£(2)) =ar, (AI9) 

and a solution exists when dV lies within the range of 
Dr- Using V^^ to denote the pseudoinverse of Vr and 
the fact that X'y(X^) = Vyixy we rewrite Eqs. (|XT7l) 

as 

L(2)=p-i(ar) (A20a) 
L^^) =2T-^dd-2V-\dT)d (A20b) 

(A20c) 

Putting this back into Eq. (jA2p we get 

£9 -:Pr^(ar),j (i?'' - d') o - d^) (A21) 
+ 2(9dT".^(i?^' -d^') (A22) 



-tr[Pr'(ar)r] 



(A23) 



This proves Eqs. ^ and 

As suggested above Vr may be singular. To ana- 
lyze the spectrum of Pr it is convenient to consider 
r — STthS"^ in Williamson form, where Fth is a (ther- 
mal) Gibbs state. In this representation. 



Vr = {S S)Vr,^{S (g, Sy 



(A24) 



and the eigenvalue equation Vr^^^{X) = \X reads 
rthATFth — ujXuj'^ = XX. Since Fth is diagonal and pro- 
portional to the identity in each 2x2 block corresponding 
to each mode, the eigenvalue equation decouples into 
independent equations, labelled by the pair cor- 
responding to the i-th and j-th modes in the row and 
column rcsp. 



UjXijLU^ = {ViVj - \)Xi^ 



(A25) 



The map ljj ■ ljj^ is an involution, thus its eigenvalues 
are ±1. Hence, the eigenvalues of Vr^^^ arc of the form 
A = ViVj ^ 1, where Vi arc the symplectic eigenvalues 
of F. In addition, matrices 1 and a; are eigenvectors of 
uj ■ corresponding to eigenvalue -1-1, and Pauli matri- 
ces ax and az correspond to eigenvalue —1. Hence l,a; 
are eigenvectors of Vr^^^ with eigenvalue ViVj — 1, and 
Oxi^z arc associated to viVj -\- 1. Thus, the spectrum of 
Ppth is given by = ViVj ^ \. In particular, one can 
see that if F has symplectic eigenvalues equal to 1 (or 
F -|- is singular), then Vr^-, is singular, and so is Vr- 
The kernel of "Drth consists of sparse matrices populated 
only in blocks (i,j) such that i^i = i^j = 1, with entries 
of the form oiij \ + ftj W, 



kcr(Dr,J 











(A26) 



0/ 



Thus Vr is nonsingular only when the symplectic eigen- 
values of F are strictly greater than 1, and if Vr is sin- 
gular, its kernel consists of matrices with with nonzero 
entries in blocks corresponding to vacuum modes in the 
Williamson decomposition of p. 



Appendix B: The quantum Fisher information 

The QFI is defined as Iq = tr[pgCg]. Using the nota- 
tion introduced in the text, a map X — > AXB^ reads 
B and the action is specified as follows 



A®B\X) ^ \AXB'^ 



(Bl) 



With this Vx = X ® X 



vijkl 



Lu has 4 upper indices 



ji-^ijjji of which the last two are 



7 



contracted with the argument matrix, 



On the other hand, supposing that T>x is nonsingular, 
[Vx%u[Vxf''''Yr, = Y,j, and 



[V-^\Y)],, = [V-^%mY'' 



and 



Tlj _,,,ikY,,,,Jl] 



= \XYX ' - ujYlo 



(B2) 



I A;/7's 



(B3) 
(B4) 



We wiU use the relation [lOf 

where the centered canonical operators R ^ R — d have been defined. With this, Iq reads 
Iq ^V^\dT),,V^\dT)kitv[pgiR' o W) o o i?')] 

+ (4[r-i5d]4r-iad]j - tr[pp i(ar)r]2?f ^(ar)y) tr[pe i?* o t?^] + ^ti-[D^^{dr)r]^ 



(B5) 

(B6a) 

(B6b) 
(B6c) 



Since ^(dT) is symmetric, we can write 

+ 2dd^T-^dd. (B7) 

Finally, notice that 

[2?r]'-'""[Ppi(9r)],, = [VroV^\dTp = [T'(ar)]'-'" 

(B8) 



I 

To conclude, notice that 2?r is self adjoint, i.e., 
tr[X'r(^)^i^] = iT[X'^VT{Y)], as a consequence of T be- 
ing symmetric. So is 2?p^. Also, the presence of 
makes V unnecessary, i.e., f p ^ oV — I'p ^, hence we can 
write 



Iq = ^tr[arPpi(ar)] + 29(i^r-i5(i 

or, in braket notation, 

Iq = \{dT\V^^\dT) + 2dd^T-^dd. 



(BIO) 



(Bll) 



where V = Vyo'D^^ is. the projector onto the support of More exphcitly, 
Vy-- Hence, 



Iq = itr[2?p 1 (9r) V{dT)] + 2dd^T-^dd (B9) 



Iq = ^(c^r|(r(g)r - uj®oj)-^\dT) + 2dd^T-^dd, (B12) 



which is the expression reported in Eq. ([T 



